Stability of ferromagnetism in the Hubbard model on the kagome lattice 
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The Hubbard model on the kagome lattice has highly degenerate ground states (the flat lowest 
band) in the corresponding single-electron problem and exhibits the so-called flat-band ferromag- 
netism in the many-electron ground states as was found by Mielke. Here we study the model 
obtained by adding extra hopping terms to the above model. The lowest single-electron band be- 
comes dispersive, and there is no band gap between the lowest band and the other band. We prove 
that, at half-filling of the lowest band, the ground states of this perturbed model remain saturated 
ferromagnetic if the lowest band is nearly flat. 

PACS numbers: 75.10.Lp,71.10.Fd 



Numerous studies have been made on the Hubbard 
model, a tight-binding model of electrons with on-site 
interactions, to understand mechanisms for ferromag- 
netism in itinerant electron systems in a simplified sit- 
uation H |, |, §. Recently, Mielke § and Tasaki § 
brought a significant breakthrough in the field by prov- 
ing that certain classes of Hubbard models have fer- 
romagnetic ground states. These models in common 
have multi single-electron bands containing dispersion- 
less bands, and are called flat-band Hubbard models. 

Although these flat-band Hubbard models shed light 
on the role of the Coulomb interaction in generating 
ferromagnetism, the models with completely flat bands 
are still singular. It is desirable to clarify whether the 
flat-band ferromagnetism is stable against perturbations 
which turn the flat bands into dispersive bands. As 
for Tasaki's version of flat-band Hubbard models, local 
stability of the ferromagnetic ground state in perturbed 
nearly- flat-band models was proved jjj . Tasaki also gave 
a concrete example of nearly-flat-band Hubbard models 
in which he could prove that the ground states are ferro- 
magnetic j|, ||] . See |l(| |ll| for related results. 

As for Mielke's version of flat-band Hubbard mod- 
els, on the other hand, there have been no rigorous re- 
sults about stability (or instability) of ferromagnetism 
in perturbed nearly- flat-band models jl2). Here we note 
that there are essential differences between Mielke's and 
Tasaki's models. Mielke's models have simple structures 
where all the lattice sites are identical, while Tasaki's 
models have two different kinds of lattice sites. Reflecting 
the lattice structures, there are no band gaps in Mielke's 
models while there are finite band gaps in Tasaki's mod- 
els. We stress that the problem of stability of ferromag- 
netism is much more subtle and difficult in Mielke's mod- 
els, where one might encounter various low energy exci- 
tation modes which arise from the gapless nature of the 
band structures. 

In this Letter , we treat the model obtained by adding 
hopping terms to the Hubbard model on the kagome 
lattice, a typical example of Mielke's models 13 . The 
added perturbation destroys flatness of the band, but 
the band structure remains gapless. We prove that our 
model has saturated ferromagnetic ground states at half- 



filling of the lowest band, provided that the lowest band 
is nearly flat. 

Definition. We first define the reference triangular lat- 
tice £ as 



£ = 



n\V\ + n 2 v 2 



ni G Z and < ni < L 
for i = 1,2 



, (1) 



where v\ = (1, 0), v 2 = (|, 4)i and L is a positive inte- 
ger. For each a G £ we define 

Cl = {x = n x f + n 2 f | m,n % € Z,|i -a\ = ±}, (2) 



C« = i x = n i Q + n 2 I ni,n 2 £ Z,\x - 
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and C a — C\ U C\ Q . Then, the kagome lattice A can 
be constructed as A = U a £cC a , where a site x G A is 
generally counted four times in different C a . (See Fig|l|.) 
Our lattice has open boundaries. Periodic lattices can 
be also treated with extra technical complication. We 
denote by <E> the state with no electrons and denote by 
c x ,a and cj. {j the annihilation and the creation operators, 
respectively, of an electron with spin a at site x in A. 
These operators satisfy the usual fermion anticommuta- 
tion relations. The number operator of an electron with 



spin a at site x is defined elS Tlx c 



= A 



The to- 



tal spin operators Stot = (^tot> ^tot? ^tot) are defined as 
S tl = \ T,xeA £<t,t=U c i,<yP^rCx,T for i = 1,2, 3, where 
p^> = [p^] a ,T=^ ,i are the Pauli matrices. We denote by 
Stot(St;ot + !) tlie eigenvalue of (S'tot) 2 - 
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FIG. 1: (a) Local lattice C a = C\ U Cl. The open and gray 
circles represent sites in C\ and C\, respectively, (b) Lat- 
tice A. 
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Let us define our Hubbard Hamiltonian on A. 
First, for each x G C\ we define fermion operator 



"x,cr ~\~ E 



yeCi;|y-x|: 



.1 C 



We also define 



a a,<T — E^ec 1 mI^j 2/] c y,o- f° r ot £ L, where coefficients 
fi[a,y] take either +1 or —1 and are chosen so that 
/j,[a,y]/j,[a,y'] — —1 whenever \y — y'\ = |. To each 
C a we associate the local Hamiltonian 

t 



'Ha — ~S 2J a a,a a a 

16C 



3 X/ X/ ^ , (Q,a;), ( T^(a,a:)^ 



(4) 



where s, t and {/ are positive parameters. Then, the 
Hubbard Hamiltonian on the whole lattice A is defined 

Remarks. It is possible to rewrite TL into the standard 

form as H = £ CT Ea^eA ^J/4, CT %,<r + E^eA ^^,T n ^,l: 
where the model parameters are given by U x — U , and 
ir,, = 2(i-s) if x = y, t xy = t+sif \x—y\ = 1/2, = -s 
if |ar - j/ 1 = V3/2, ^ = s if \x - y\ = 1 and x,y G C£ 
with some a, and i xy = otherwise, except for the sites 
close to the boundary. The single-electron dispersion re- 
lations (calculated in the model with periodic boundary 
conditions) are given by E\{k) — — 2s(3 — e(k)), ^(fe) = 

i(3 - V 3 + 2e O))> and ^(fc) = t(3+ y/3 + 2e(k)) 
with e(fc) = cosfci + cosfc2 + cos(fci — /c 2 ), where 
fc = fci^J + ^2^2 is the wave vector expanded in 
terms of reciprocal- lattice vectors z/* = (l,— and 

i/* = (0,^). Note that £i(0) = £ 2 (0) = 0, which 
means that there is no gap between the lowest and the 
second lowest bands. 

One readily finds that {a Q „, &L ^ a } = for any 

a,fl G £ and a; G C|. This implies that {a^ a ^ } ae c 
spans the space corresponding to the lowest band. 

If we set s = 0, our model has highly degenerate single- 
electron ground states, and becomes essentially the flat- 
band model of Mielke's (although there is a difference in 
boundary conditions). In this case, the model exhibits 
flat-band ferromagnetism for all positive values of U. In 
the model with s > 0, the situation is quite different 
because double occupancies of lower energy states, which 
destroy the ferromagnetic order, may reduce the total 
energy of the system. It is indeed easy to prove that 
the ground states of our model has 5 t ot = (or |) for 
[7 = 0, and cannot exhibit saturated ferromagnetism for 
sufficiently small U. (See, for example, Sect. 3.3 of Q). 
The following theorem establishes that the ferromagnetic 
ground states are stable for sufficiently large t and U 
when the electron number is \C\. 

Main theorem. Consider the Hubbard model defined 
as above with the electron number \C\. Then, there ex- 
ist critical values (t/s) c and (U/s) c , independent of the 
lattice size, such that, if both t/s > (t/s) c and U/s > 
(U/s) c are satisfied, the ground states of the model have 
Stot = Furthermore the ground state is unique up 

to the degeneracy due to the rotational symmetry. 



In Tasaki's models, the stability of ferromagnetism 
may be, at least at a heuristic level, understood as a con- 
sequence of the band gap separating the lowest nearly 
flat band from other bands. The band gap enforces the 
electrons to occupy the lowest band while the interaction 
rules out double occupancies of sites. Then the situation 
is almost as in the flat-band models, and the systems ex- 
hibit ferromagnetism. To Mielke's models, which have 
no band gaps, the above argument does not apply, and 
the origin of the stability of ferromagnetism seems more 
subtle. Nevertheless, our proof is based on essentially the 
same philosophy as that of Tasaki's proof in || . Namely, 
we first establish ferromagnetism in a local model de- 
scribed by TL a , and then show that these local ferromag- 
netism can be "connected" , which results in macroscopic 
ferromagnetism in the whole system. The results of the 
analysis of TL a are summarized in the following lemma. 

Lemma. If t/s and U/s are sufficiently large, the min- 
imum eigenvalue of 7i a is —6s and any eigenstate $ be- 
longing to this eigenvalue is written as 



$ = a 1 t $t + a 1 <E> 

ct,T I a, I * 



(5) 



with appropriate states $j and $| . Furthermore it satis- 
fies 



c s ,|C X)T $ = 



(6) 



for all x £ C a • 

Proof of Lemma. Since all the local Hamiltonians are 
the translated copies of 7io, it suffices to prove the lemma 
for a = 0. From now on, for convenience, we identify Cq 
and Cq with {0,2,..., 10} and {1,3,..., 11}, respectively 
(we first label (4,0) as 0, then label the rest sites as 
1, . . . , 11 in the clockwise order). 

We start by solving a single-electron problem for Ho- 
Let I = {0, ±-|, ±Tp 7r}. Then the eigenvalues are given 
by 



Slip) 



ifpeAM; 
—6s if p = 7T , 



(7) 



and £2(p) = |(3 + 2cosp) with p G I. The eigenstate 
corresponding to £i(p) is expressed as dp a $>o with 



dp. a 



1 



E^ 



— ^ [C2l,a ~ C 2 l-l,a - C 2 l + l.a) 

a/6(3 + 2cosp) £5 

(8) 

(where c_i, CT is regarded as C\\ Note that the set 
{d+ CT $o}peJ is orthonormal since {d PjCr , d p , a } = S p y. 

We consider many-electron problem for 7io , first in the 
limit t,U —^00. Let $ be a state on Co which has a finite 
energy in this limit. Since all e-i[p) are infinite in the limit 
t — > 00, $ must be expanded as 



(9) 



IrJlCl 
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with complex coefficients g{I-\, where 

$(j T; / i )= n< T n 4,^° • ( io ) 

pel-, p'eh 

Here, and throughout the present Letter, the products 
are ordered in such a way that d} p ^(resp. <ip i) is always 
on the left of d p , ^(resp. dL , ) if p < p'. Since the on-site 

interaction n x ^n x ^ — c x *c x ^c x ^c x ^ is positive semidef- 
inite, the state $ in the form of (|J) must further satisfy 

E g(I r ,I l )c x , i c xA $(I r ,I i ) = (11) 

7 T ,7 iC 7 

for any x G Co in order to have finite energy in 
the limit U ->• oo. From (0), (|) and (0) one 
finds that the expectation value of 77o for the state 
$ is = ($,-H $)/(*,$) = -6s + 6sF||$|r 2 , 

with F = £/ T ,/ iC /\ W (\g(A;h)\ 2 - l3(If;/f)| 2 ) and 

ll^ll 2 = (E/ T ,7 iC 7l5UT;^)l 2 )[El: where coefficients g 
should satisfy the condition (pi]). Here, and in what fol- 
lows, we abbreviate I a U {p} as 7£ for pel. In the 
following, we will show F > 0. This implies -E<j> > —6s 
since s > 0. 

To prove 7 1 > 0, we first derive conditions on g imposed 
by eq. ©. If we denote (y4 p) )* = {c a , CT , g^J, the left- 
hand-side of eq. ( |TT| ) becomes 

E sUt;^) E E (-i) |/Thl s? T S^ 

7 T ,7 i c7;|7 T |>i,|7 1 |>i pehp'eh 

= E (^ p) n^ p,) r 

p,p'G7 

x E E (-i) |/T| s? T S?>(if;/f>(/T;^) 

/ T c/\{p} ^cz\{p'} 

= E E (^ p) )*(^ p ' ) )^( / f; / r>( / T^i) ! (12) 

ixJiCi p,p'ei 

where Sj , which corresponds to a sign factor com- 
ing from exchange of the fermion operators, equals 1 
if S P '6/„y<pl i s even an d ~1 otherwise. In the fi- 
nal expression of (|l2[), wc introduced subsidiary co- 
efficients g defined as g(I^;I^) = if p G 7j or 

p> G 7 4 and 5(/f;/f ) = (-1)^S? T S^T*; if) oth- 
erwise. Therefore, c x ^c x ^<& = holds if and only if 

E P ^ e /(^ P) )*(^ P ' ) )*5(/ T P ;/f') = for any 7 T ,7 4 C i. 
Taking the sum of this equation over x G Cg, we 

find that Epgj (3+2cosp) g( / fi / r P ) = and simi " 
larly taking the sum over x G Cg, we find that 

E pe j ( 3+2^ P p) g( J f ; j r P ) = ( whcre we identined 
with tt). By eliminating g(l9;I?) from these two equa- 



tions, we obtain 

-^(/ T ¥ ;^ ¥ )-^f ¥ ;^).(i3) 

Our analysis below relies heavily on this condition. 

For a subset 7j of 7, wc define 7j_ = {— p | p G TjJ 
and denote by iV(7|-;7|_) the number of elements in 
7 T n 7 X n (7\{0,tt}). Condition @ relates §(7 T ;7 i ) with 
7 T ,7 X such that A?'(7 T ;7 i ) = r and g{I\;I{) with 7f,7j 
such that AT(7|;7{) = r + 1. This motivates us to decom- 
pose F as 7 1 = 7 1 ' + ^ =Q F r , where 

F r = Wr,h)\ 2 - E Wl^i)W 

I u hCl\{7r}; 7 T ,Z 1 c7\{ 7 r}; 
JV(/ T ;Z 1 )=r+l JV^^r 

(14) 

= E ls(^i)l 2 - (is) 

J T ,7 i CJ\{7r};JV(/ T ;Z i )=0 

Since the term F' is apparently non-negative, F > is 
implied by F r > for r = 0, . . . , 4. 

We shall prove that F r > by using ( |l3| ) . For a pair of 
i? and 7^ such that N(II; 17) = r, the number of non- 
zero g in the right hand side of (jl^) is, by the definition, 
at most 4 — r, and thus for such a pair we have |]l6| 

l5(/f;^)l 2 <^( 4 -^) E Iff^fT- (16) 

pG7\{0,7r} 

Then, we find that 

E ls(/W)l 2 = E W;^)! 2 

Ar(/ T ;/;)=r WC/^J^)^ 

^ £(*-r) E E i^ T vr p )i 2 

/ T ,7 i C/\{7r};pe/\{0, 7r } 
W(I T ;ii)=r 

= ^( 4 -^ + 1 ) E l.9(^ T ;^)l 2 

7 T ,7 i c7\M; 
iV(7 T ;Z 1 )=r+l 

< | E I5(a;^i)I 2 - (i7) 

-fT,-fiC7\{7r};AT(7 T ;7 i )=r+l 

To get the third line, we have used the fact that, for 2j 
and 7| such that 7V(7| ; 7^) = r+1, there are r+1 elements 
p in 7\{0, 7r} for which we can find suitable subsets 7^ and 
I[ such that {p} U II = 7 T and {-p} U 7| = 7|. To obtain 
the final inequality, we have used (4 — r)(r + 1) < 6 for 
< r < 4. By using (jl^) we obtain 

F r >^ E |.g(7 T ;7 i )| 2 >0. (18) 

7,,7 i c7\{7r};Af(7 T ;7 i )=r+l 



4 



We therefore conclude that F > 0. The above analysis 
also shows that the equality F = holds only when F' 
and F r are vanishing, i.e., g(If ; ij) = for any pair of If 
and 1 1 such that 7r € If PI Ij or 7r ^ 7| U 7| . 

In other words we have shown that E$, > —6s for any 
$ and that any $ attaining the minimum expectation 
value —6s is written as 



4> 



E 



(19) 



If , 1 1 C 1 1 IT G 7 ? U 1 1 , 7T £ If n 7; 



and further satisfies the finite energy condition (^lj) . One 
finds that such minimizing $ indeed exists by testing 

4, T $ o or n 



c£ T $o- By construction such $ is an 



pel p,T 

eigenstate of Ho as well as — s aj a ao i(J . Since it is 
known to be the lowest energy state of Ho in the limit 
t, U — > oo, the continuity of energy implies that such <f> is 
the lowest energy state of Ho for sufficiently large t/s and 
U/s. It is also easy to check that such $ has the prop- 



x/6 



ao,ff-) 



erties stated in Lemma. (Note that 
This completes the proof of Lemma. ■ 

Proof of Theorem. We assume that the values of t/s 
and U/s are large enough for the statement in Lemma to 
hold. We note that how large t/s and U/s should be is 
independent of the size of A, because Lemma is concerned 
with the local Hamiltonian. 

From Lemma wc find that the eigenvalue of H 
is bounded below by — 6s|£|, while, by taking 

$ f = n 



<E>o as a variational state, we find that 
-6s|£| is an upper bound on the ground state energy. 
Therefore, the ground state energy is — 6s|£|, and $f and 
its SU(2) rotations are among the corresponding eigen- 
states. It is apparent that these states have S to t = |£|/2- 



The remaining task is to prove the uniqueness. Let 
$G be an arbitrary ground state of H. Lemma implies 
that the ground state energy is attained if and only if 
H a $G = — 6s$ G for all a e C. Thus $ G must satisfy 
the conditions stated in Lemma. 

The condition (j^) implies that <E>g is expressed as |l7j 



*G = £V(M)II U.*°: 



(20) 



M 



where {a} is a shorthand for a spin configuration 
{<7 a } , 6 £, the summation is over a a =|, J, for all a e C, 
and ip({a}) is a complex coefficient. 

Let us impose the condition (^J) on <I>g in the form 
of (20). Let j3 and 7 be nearest neighbour points in C, 
i.e., 1/9 — 71 = 1, and let m(/3, 7) € A be the site located 
at the mid-point between f3 and 7. It is easy to see that 



i(/3,7), 



r } is nonvanishing if a = [3 or 7, and is 



vanishing otherwise. Then, it follows from the condition 
c m (/3, 7 ),|Cm(/3, 7 ),T$G = that v?({cr}) = ¥>({t}) for any 
pair of spin configurations {a} and {r} satisfying that 
o~/3 = Try, cr 7 = 7/3, and cr Q = T a for a ^ /Jj. Examin- 
ing c m ( / 3 :7 ) i |C m ( | a )7 ) ! i-$G = for all the pairs of nearest 
neighbour points in C, we find that ip({a}) — ip({r}) 
whenever Y] a o~ a — ^2 a T a . Therefore <!>g is written as 
o i PM(S tot ) M ${ , where </?m are new coeffi- 
cients and the spin lowering operator 5 t ~ t is defined by 



^tot — Yl 



This completes the proof. 
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